Comment on " Structure of a Quantized Vortex near the BCS-BEC 

Crossover in an Atomic Fermi Gas" 
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Machida and Koyama presented a study of the quan- 
tized vortex core structure near the BCS-BEC crossover 
regime [1]. A very similar analysis was performed by 
Kawaguchi and Ohmi [2]. The conclusions are similar, 
based on essentially the same theoretical approach due 
to Timmermans et ai, [3]. Initially this theoretical ap- 
proach was believed (incorrectly) to handle in a satis- 
factory manner the case of the large scattering length a, 
when n\a\ 3 > 1, where n is the atom number density. In 
this approach one introduces a boson degree of freedom, 
associated with two atoms forming a boson molecule in 
the closed channel. The practitioners of this approach, 
conclude typically that in the BCS-BEC crossover regime 
there is a significant, even dominant, probability of the 
atoms to be in the closed channel. If that would be the 
case, then an atomic Fermi gas in this regime would be- 
have essentially as a Bose system. The probability of 
being in the closed channel is negligible [4] and a re- 
cent direct measurement of this quantity confirms this 
unequivocally [5]. The authors of Refs. [1,2] seem to ar- 
rive at qualitatively similar conclusions to those obtained 
by the author and Y.Yu in Refs. [6] , that in the BCS-BEC 
crossover regime a vortex in an atomic Fermi gas shows 
an unexpected [7] marked density depletion. However, 
these authors argue, incorrectly, that this density deple- 
tion is due (mostly) to the significant, even dominant, 
occupation probability of the closed boson channel. 

On one hand, the two channel approach [3] suffers 
from several deficiencies. While the physics is clearly de- 
termined by a single dimcnsionlcss parameter only, the 
typical one used being l/kpa, where n — kp/3% , the 
model of Refs. [1-3] is overdetermined (U - the atom- 
atom "bare interaction," g - the closed channel boson to 
two atom coupling, v - the detuning of the closed channel, 
an ill defined energy cutoff, etc.). Moreover, the typical 
usage of this model is within meanfield (with some fluc- 
tuations included sometimes). It is well known that the 
corrections to the meanfield are controlled by the param- 
eter kp\a\ ^> 1. Moreover, a certain type of fluctuations 
(which are routinely ignored in these treatments) lead to 


a strong reduction of the pairing field both in the weak 
[9] and strong coupling limits [8]. The calculations of 
Refs. [1,2] neglect the role of the attractive meanfield 
too, which thus disfavors a density depletion. More to 
the point, experiment [5] shows unequivocally that in the 
BCS-BEC crossover the boson component contributes 
« 3 x 10~ 6 • • • 2 x 10~ 4 , as opposed to the theoretical 
predictions of Refs. [1,2], namely « 0.4- ■ ■ 1. Clearly, 
such an insignificant (as observed) fraction of (compos- 
ite) bosons cannot influence the vortex core structure. 

On the other hand, the calculation of Refs. [6] are 
based on a theoretically consistent extension of the DFT 
[10] to fermionic superfluid systems [11] and on fully non- 
perturbative calculations of the homogeneous state [8], 
defined by one parameter only, namely l/kpa. These 
results show, that in spite of the quenching of the pair- 
ing gap due to fluctuations, and in the absence of any 
boson contribution (which would otherwise favor a den- 
sity depletion), there is a significant density depletion at 
the vortex core, qualitatively consistent with the recent 
experimental observations [12]. 
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